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I discuss topics related to four-quark states of tlie 'molecular quarkonium' type, i.e. resonances that could be 
considered as (dominantly) made from a heavy meson and antimeson. Of the so far observed resonances such 
picture is very likely applicable to the state X(3872), and I also discuss its possible relevance to the peak near 
the D* D* threshold in e+e~ annihilation. 
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1. Introduction 

It has been pointed out long time ago ,1] that 
there should exist bound and/or resonant states of 
heavy hadrons, dubbed 'molecular' and arising due to 
the strong force between the light constituents of the 
heavy hadrons. Indeed, considering the interaction 
between a heavy meson, consisting of a heavy quark 
Q and a light antiquark g, and an antimeson {Qq), it 
is clear that there is a normal strong force between 
such objects arising from an exchange of a light qq 
pair, which can be thought of as an exchange of light 
mesons. The range and the strength of the resulting 
strong interaction does not depend on the mass mg of 
the heavy quark Q in the limit mg — > oo. If the inter- 
action corresponds to an attraction in some states of 
the heavy meson- antimeson pair, it should inevitably 
give rise to bound states in those channels at suffi- 
ciently large toq. However, given the lack of detailed 
knowledge of the strong interaction between heavy 
hadrons, it has not been clear, whether the charmed 
quark is heavy enough for the reasoning based on the 
heavy quark limit to be applicable, and thus whether 
resonances of such type with hidden charm do exist. 
Another difficulty of identifying states with hidden 
charm as molecular, or more generically as four-quark 
systems, is that in most models of light meson ex- 
change between charmed mesons 0,0,13 the interac- 
tion with the strongest attraction arises in channels 
with zero isospin and with non-exotic quantum num- 
bers J^*^, which can be also realized as a conventional 
cc charmonium. For this reason it could be a typi- 
cal situation that some dominantly four-quark states 
may be indistinguishable from states of charmonium, 
unless some other properties would favor a molecular 
interpretation. One such interpretation was argued in 
the past lUi for the resonance in the e+e" annihilation 
at the threshold of the D*D* production, then called 
V'(4028), and presently referred to Q as ip{4040). It 
was noted that the cross section for production of the 
D*D* pairs accounted for more than one half of the 
total production of charmed meson pairs in spite of 
the very strong suppression by the P wave kinemati- 
cal factor at energy close to the threshold. Thus it 
was suggested that there is a resonance very close to 
the threshold actually made from D*D* meson pairs. 



The idea of molecular states was revived in a 
modern context by the observation |^ in the Belle 
experiment of the resonance X(3872) decaying into 
7r+7r~J/^/;, and produced in the decays B ^ X K of 
the B mesons. The resonance is quite narrow, and 
its width is not yet resolved {Tx < 2.3 MeV at 90% 
CL), and the updated value of its mass |7J mx = 
(3871. 2±0. 6) MeV coincides within the errors with the 
mass of the D°D*° meson pair: Af(L>°) + M{D*") = 
3871.2 ± 1 MeV. The extreme closeness of the newly 
observed resonance to the meson pair threshold is very 
suggestive of an interpretation [EM 113)113 of this res- 
onance as a 'molecule' dominantly consisting of neu- 
tral charmed meson pairs: D°D*° ± D*°D°. 



2. X(3872) 

The discovery and observation 0, 0, 0| 
of X(3872) through the decay X tt+tt-J/^P 
was followed by finding also the decays X 
7r+7r^7r° J/?/; and X iJ/"^- Clearly, the co- 
existence of the latter decays with the discovery 



mode X —> tt+tt^ J/i/) implies an isospin breaking in 
X(3872). Indeed the two-pion state and the three- 
pion one have opposite G parity, and one or the other 
would have to be forbidden if the isospin symmetry 
were applicable in the decays of X(3872). Further- 
more, the existence of the radiative decay X "fJ/ip 
requires X to have positive C parity, and the two- 
pion system emitted in the decay X —>■ tt^tt^J/^ 
has to be in a C-odd state, which corresponds to the 
isospin of this system equal to one: /tttt = 1- Thus 
the X(3872) resonance is required to have a substan- 
tial 1=1 component, which is impossible for a pure 
cc charmonium. In other words, the resonance has 
to contain light uu and/or dd quark pairs in addition 
to hidden charm, which thus qualifies it as a four- 
quark state. Moreover, the decays X — ^ tt+tt^ JM 
and X 7r+7r~7r'' J/i/) have comparable rates p!q |. 
which implies that the X(3872) resonance does not 
have a definite isospin, and is mixture of / = 1 and 
/ = states. In addition, a subsequent study 
of the angular correlations and the dipion mass spec- 
trum in the decay X —^ n'^'K^J/ip strongly favors 
the spin-parity assignment J = 1+ for the X(3872). 
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Combined with the extreme proximity of its mass to 
M(D°) + M{D*°) the known properties of X(3872) 
most strongly suggest that a sizeable part of the wave 
function of this resonance is contributed by the C-even 
S'-wave state of charmed meson pair D^D*^ + D*^D^, 
which thus has the quantum numbers J^*^ = 

The binding energy w of the neutral meson pair in 
X(3872) is small: likely w<l MeV. In this scale the 
mass gap d between the corresponding threshold for 
the pair of charged mesons and that for the neutral 
ones 6 = M{D+) + M{D*+) - M{D") - M{D*") « 
8 MeV is large enough to suppress the weight of 
the state D^D*~ + D~ D*~ in the wave function of 
X(3872), thus explaining the unusual isotopic proper- 
ties of the resonance. 

Generally, the wave function of the resonance 
X(3872) can be written in terms of a Fock decom- 
position: 

\X) = aQ\D^D^* +D^*D^) 

+ ai \D+D-* + D+*D-) + lothcr) , (1) 

where first two terms correspond to the 'molecular' 
states of charmed meson pairs, and the "other" states 
in the sum describe the admixture of non-molecular 
states, such as e.g. a ^Pi charmonium cc pair. The 
latter states are localized at shorter distances within 
the range of the strong interaction, while the 'molec- 
ular' part refers to the motion of the charmed mesons 
at (mostly) the distances beyond the range of strong 
interaction [Tol[l7l |. At the latter distances the coor- 
dinate wave function (f>(r) for each pair is that of a 
free motion: 



distances the wave function of X{3872) is dominated 
by / = 0, so that the 'core' of the wave function of 
X is an isotopic scalar. Then the wave functions of 
the D+D*- + D-D*+ and D°D*° + D*°D° should 
coincide at short distances, which allows to estimate 
the relative statistical weight of these two states in the 
X(3872) resonance (i.e. the ratio of the corresponding 
coefficients in Eq. QJ: 



A = 



|ao| 



(4) 



Strictly speaking, the wave functions in Eq. [21 and 
Eq. 131 cannot be applied at short distances in the re- 
gion of strong interaction, where the mesons overlap 
with each other and cannot be considered as indi- 
vidual particles. In order to take into account this 
behavior an 'ultraviolet' cutoff should be introduced. 
One widely used method for introducing such cutoff 
is to consider the meson wave functions only down 
to a finite distance tq, at which distance the bound- 
ary condition of the state being that with / = is 
imposed (a discussion of a similar situation can be 
found e.g. in Ref . [T^ ) . An alternative, somewhat 
more gradual cutoff, described by parameter A, can 
be introduced^3 by subtracting from the wave func- 
tions Q and ||3Jl an expression ce~^^ /r. One can also 
readily see, that an introduction of any such cutoff 
eliminates relatively more of the charged meson wave 
function than of the neutral one, thus reducing the 
estimate of the relative statistical weight as compared 
to that in Eq. 0] so that Eq. 0] gives in fact the upper 
bound for the ratio. 



(Pn(r) = J- 

V ZTT r 



(2) 



for the meson combination {D°D°* +D"*D°)/ ^2 and 

(3) 



/«~exp(-Kcr) 
Mr) = \ 7y 

V ZTT r 



for {D+D-* + D+*Z)-)/^/2. The virtual momenta 
K„ and Kc are related to the binding energy w and 
the reduced mass of the two-meson system w 
970 MeV as k„ = ^2 rurW « 44 MeV ^Jw/MeW and 
Kc = ^/2rn7W+w) « 125 MeV. 

The dominance of the D°D*° + D*°D° at long dis- 
tances translates into a substantial isospin violation 
in the processes determined by the 'peripheral' dy- 
namics, examples of which are apparently the ob- 
served decays X Tr~^Tr~ J/ip and X tt+tt^tt^ J/?/;. 
It is quite likely however that this isospin-breaking 
behavior is only a result of the 'accidentally' large 
mass difference 5 « 8 MeV between D^D*^ and 
£)*Ojj*o Therefore it is natural to expect that at 
shorter distances within the range of the strong inter- 
action the isospin symmetry is restored and at those 



3. Peripheral Decays of X(3872) : 

X D^D^n^ and X ^ DO^O' 



7 



The peripheral D°D*^ + D*°D° component with 
the quasi-free mesons should give rise to the decays 
X D°D°Tr° and X D°D°j due to the underly- 
ing decays of the D*'^ and mesons. The widths of 
the corresponding underlying decays can be deduced 
from the available data @: E^ = r{D*° D^tt") = 
43 ± lOKeV and = T{D*° D°-f) = 26 ± 6KeV. 
In the C-even state X(3872) the amplitudes of the 
decays D*° D°Tr° and D*° D°tt° interfere con- 
structively and the result for the rate of the decay of 
X can be written as [l^] 

TiX^D'^D\°) = 2\ao\'T^[A{w)+7^B{w)] , (5) 

with A{u)) describing the non coherent contribution of 
the decays from D*^ and D*^ and B{w) being the in- 
terference term. The dependence of these coefficients 
on the binding energy w is shown in Figure ^ As 
can be seen from the plot the interference is quite 
important at w as small as 0.1 MeV and its relative 
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importance increases with w, since at stronger binding 
the distance between the mesons in X(3872) becomes 
shorter. The overall fall-off of both A and B at larger 
w is the result of the decrease of the phase space for 
the decay. 



dr/dffl 




Figure 1: The dependence of the coefficients A and B in 
Eq.|3on the binding energy w. 

The amplitudes of the radiative transitions D*^ 
D^^ and D*'^ D^^ interfere destructively in the 
decay X(3872) D°D°j. ^ The expression for the 
rate has the form [TlflTll: 



arctan - 



(6) 

where the photon energy lo can be approximated by 



that in the free meson decay D* 



137 MeV, due to the fact that the weak binding only 
slightly spreads the photon spectrum in the decay 
X D^D^'y near ujq, as can be seen from the plots 
shown in Figure |21 The shape of the spectrum is de- 
termined by the wave function of the meson pair in 
the X(3872) and thus a measurement of this shape 
is a direct probe of the motion of the mesons in the 
'molecule'. 

It can be noted that the other kinematically al- 
lowed radiative decay channel X D^D^j should 
be quite suppressed [13 ■ The suppression of the di- 
rect peripheral contribution due to the transitions 
7 arises from the discussed small over- 
all statistical weight of the D+D^* + D+*D^ state 
in the wave function and from a strong destruc- 
tive interference between the amplitudes of the pro- 
cesses D*~^ — > Z?+7 and D*~ D~j due to a sig- 
nificantly shorter spatial separation of the charged 



^Such effect of destructiv e interfe rence in the photon emis- 
sion is known is subradiance 113,123. 
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Figure 2: The photon energy spectrum in the peripheral 
decay X(3872) D°D°'y, at w = 0.3 MeV (dashed) and 
at u) = 1 MeV (solid) 



mesons in X(3872). The contribution due to the final- 
state rescattering X — > D^D^^ D~^D~j is also 
very small due to the behavior of the scatter- 
ing D^D'^ D^D~ very near the threshold below the 
■0(3770) resonance. The most efficient way of produc- 
ing the D^D^-f final state in the decay of X(3872) ap- 
pears to be due to the transition X(3872) i/i (3770)7 
from the 'core' of the X resonance. Notice how- 
ever that such transitions result in the photon energy 
u! Ri 100 MeV i.e. well below the peak from the pe- 
ripheral process (cf. Figure 



4. The Core of X(3872) 
4.1. Spin Selection Rule 

The structure of the 'core' of the X(3872) resonance 
represented by the "other" states in the wave func- 
tion JQ) is largely unknown and in discussing this part 
we have to resort to general considerations and quite 
approximate estimates. One such general property 
based on the heavy quark limit for charmed quarks is 
the spin-selection rule for the cc pair in X(3872) j23 |. 
Namely, unlike a generic four quark state, the S wave 
C-even (D^D*" + D°D*^) state is unique in the sense 
that the total spin of the heavy cc pair is fixed: 
Sec — 1- Indeed, for the S wave there is no orbital 
motion, and the total spin of the state is given by 
the vector sum of the spins of the cc and uu quark 
pairs: J — Sec + Su-a- However the states of the type 
Sec = (X) Suu ~ 1 and Scs — I Sun — have nega- 
tive C parity: C=-l. Only the state Scs — Sua — 1 
results in C=-t-l. The total spin of the light quark pair 
is not 'traceable', i.e. can be flipped in the mixing with 
the "core" states. The heavy quark spin however is 
conserved in the limit mg oo. Therefore one arrives 
at the spin selection rule (valid up to 0{Aqcd /f^c))'- 
the states with Sec = 1 dominate in the Fock sum in 
Eq.HJfor X(3872). 
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As one application of the spin selection rule one 
can consider hadronic transitions from X(3872) to the 
charmonium states. The spin selection rule implies 
that the transitions to the spin-singlet charmonium 
levels such as X(3872) nnijc should be suppressed 
while the transitions to spin-triplet states should be 
allowed. The latter include the observed decays X 
7T~^Tr~ J/t/j, X 7r"'"7r~7r° J/?/) and X jJ/ip, and 
also the yet unobserved decays to the spin-triplet Xcj 
levels: X tt'^XcJ- The conservation of the spin of 
the cc pair implies the statistical behavior of the decay 
rates: 



B{B+ XK+). The current data 2|| on these de 
cays 

B{B° XK^) 



V{X ^ ^°xcj) « (2 J + 1) 



(7) 



and an approximate estimate of the absolute rate [22l| 
suggests that the width of the decay X — > tt^XcI can 
be as large as about 0.3 r(X -k^ ir^ J / tj:) . 

4.2. Production of X(3872) 

The core of the wave function of the X(3872) reso- 
nance is located at short distances. Thus it is natu- 
ral to expect that the production of the resonance in 
hard processes proceeds through the production of the 
core. Furthermore, the quantum numbers J^'~^ = 
can be realized in a '^Pi cc state, so that one can 
expect that the production properties of the X reso- 
nance should be similar to those of charmonium with 
the rates scaled by the statistical weight of the core 
jocp. The hard processes where X(3872) has been 
observed so far are the decays B XK and the pro- 
duction in the proton-antiproton annihilation. The 
analysis by CDF suggests that the production of 
the X(3872) in the pp collisions is indeed similar to 
that of the ip' charmonium resonance. The similarity 
between the production of X(3872) and of the pure 
charmonium states allows to approximately estimate 
the weight factor of the core |acp. Indeed, it is an 
experimental fact that the known charmonium states 
are produced in B decays in association with a sin- 
gle Kaon with approximately the same rate (within 
a factor of two). One might expect then that the 
core states of the X(3872) are produced in similar de- 
cays at approximately the same rate, so that the only 
suppression factor in the observed Q rate of decays 
B ^ X K \s that of the probability weight |acp, i.e. 
few percent. 

Furthermore, as discussed, the core part of the 
wave function ^ is expected to be dominantly an 
isotopic scalar, thus within such production mecha- 
nism one should expect that the decays B^ X K'^ 
and B^ X should have approximately equal 
rates 22j. This conclusion is significantly different 
from the calculation |23| of the ratio of these decay 
rates, based on the mechanism where the resonance 
is produced through its peripheral meson component. 
The latter mechanism predicts B{B^ X K°) < 



B{B+ ^ XK+) 



0.50 ±0.30 ±0.05 



(8) 



tend to favor the core production mechanism, al- 
though an improvement of the experimental accuracy 
would be certainly helpful. 



5. Models 

The finding of X(3872) and the realization that this 
is likely a dominantly molecular-type state were quite 
unexpected, and the problem of understanding the 
dynamics of such states and predicting other possi- 
ble instances of 'molecules' still remains. Although 
early papers P, 0) 01 have provided the general idea 
of molecular hadronic systems, a detailed description 
of the strong interaction involving heavy and light 
quarks is still to be developed. In this section I briefly 
review some models which have been suggested after 
the discovery of the X(3872). 

The model suggested by Swanson involves a 
calculation of the dynamics of coupled channels in- 
cluding the DD* , DD* pairs and also the pairs pJ/tp 
and ajJ/ip, assuming a certain specific form of the 
Hamiltonian for rearrangement of quarks between the 
hadrons. The model produced a successful prediction 
of the decay X — > ujJ/4> {X 7r"'"7r~7r° J/V') with the 
rate comparable to that of the discovery decay mode 
X — > pJ/ip {X TT+vr^ J/i/)). The model provides a 
detailed description of the decays of X(3872) to these 
channels as well as to DDtt, DD^, and Tr'^jJ/tp. I am 
not aware of predictions of this model for other res- 
onances of this type with hidden charm, however the 
model predicts a 1"'""'' resonance in the BB* system at 
10562 MeV and a 0"+ state at 10545 MeV. 

In the scheme of Maiani et.al. |2^ the resonance 
X(3872) belongs to a large class of four-quark states 
for which the building blocks are color-antisymmetric 
diquarks and antidiquarks. In this model this res- 
onance is identified as a 1++ state of the structure 
Xu = [cu]s=i X [cu]s=o + Ms=o X [cu]s=i. The 
model predicts also Xd = [cd][cd] resonance with the 
mass splitting M{Xd) - M(X(„)) « 8 ± 3MeV, and 
also charged partners X^. None of the predicted reso- 
nances with masses very close to that of X(3872) were 
found so far. It can be mentioned in connection with 
a possible resonance of the X^ type, that even if it ex- 
ists, it may be quite difficult to observe it, since it is 
likely to be broad due to the allowed decay into pairs 
of neutral mesons: D°D*° and D*°D. This remark 
however does not apply to the charged states X^. 

Karliner and Lipkin 27] have recently considered a 
nonrelativistic quark model with pairwise interactions 
between the constituent quarks and antiquarks de- 
scribed by an oscillator potential and proportional to 
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Figure 3: The interaction through an isovector exchange between the neutral mesons has opposite sign compared to 
that between a neutral and a charged meson 



the gluon-exchange-type color correlation within each 
pair. The Hamiltonian in the model has the form 

with Ac being the color A matrices. The authors find 
that the model predicts four-quark bound states of the 
bcqq type, but not of the ccqq one. They suggest that 
the resonance X(3872) can still be explained within 
this model, once the spin-dependent forces are taken 
into account. 

The scheme considered by Braaten and 
Kusunoki p3 | phenomenologically treats the near- 
threshold dynamics of the _D°Z)*° + DD*'^ pair as 
being exclusively and universally determined by 
the pole in the scattering amplitude, corresponding 
to X(3872), and conversely, the properties of the 
resonance as being totally determined by the meson 
pairs. This approach corresponds to an 'extreme' 
molecular treatment of the resonance. In other words, 
it assumes the limit, where any other states in the 
Fock decomposition in Eq. ^ are neglected in com- 
parison with the first term, whose coordinate wave 
function is given by Eq. [3 Clearly, such approach 
would miss the properties of X(3872) determined by 
its 'core', in particular, as discussed above, it appears 
that the core dominates the production of X in hard 
processes. 

As a general remark, it can be mentioned that exis- 
tence of the X(3872) as a dominantly D"D*° + DD*° 
shallow bound state, does not necessarily imply ex- 
istence of nearby resonances, in particular of the 
charged partners X^. Indeed, if a strong force due 
to exchange of isovector light-meson states (p, tt, . . .) 
contributes to the binding in the D'^D*° + DD*° chan- 
nel, the sign of such force would be reversed in the 
charged states, as can be readily seen from the Fig- 
ure 13 the relative sign of the two vertices for the 
exchanged isovector state is different between the two 
shown systems. 

6. D*D* Threshold in e+e Annihilation 

The very strong enhancement of the e+e~ produc- 
tion of the D*D* meson pairs very near the threshold 



is known for a long time, and it served as motivation 
for the suggestion 4] that this is an effect of a 'molec- 
ular' type P-wave J^'^' = 1 state of D*D*. In- 
deed according to the new CLEO-c measurements j2^ 
(which also agree with the old data 5]) the cross sec- 
tion for e+e~ D*D* at Ec.m. — 4030 correspond in 
units of i? to Rjjtjj, w 0.75 even though the velocity v 
of the D* mesons at this energy is u ~ 0.1 and value of 
Ro'D- is proportional to the P wave factor k, 10"'^. 
It appears quite unlikely that such a strong enhance- 
ment of the production is possible without a singular- 
ity near the D*D* threshold. Furthermore, the fea- 
tures in the cross section of the e^e~ annihilation into 
channels with pairs of lighter mesons: DD, D*D and 
DgDs around the same energy, although present, are 
still of moderate amplitude '29] . This behavior points 
to an existence of a near threshold resonance, which 
is dominantly coupled to the D*D* channel and has 
only moderate coupling to pairs of lighter mesons. 

It should be mentioned that the strong p'^ growth 
of the cross section for production of the D*D* pairs 
above the threshold leads to the behavior where the 
local maximum of either the total cross section of the 
e+e~ annihilation (located close to 4040 MeV, hence 
the notation -0(4040) or the cross section for the 
D*D* production, generally does not coincide with 
the actual position of the resonance. This position 
should be below the maximum of the D* D* produc- 
tion cross section, and can be either below, or above 
the threshold, or even between the thresholds for the 
neutral ]j*of)*o g^j^^j charged *D*^ vector mesons, 
which thresholds are split by about 6.6 MeV. 

The possible existence of a resonance dominantly 
coupled to D*D* generally leads to interesting fea- 
tures in the cross section for production of lighter 
meson pairs due to channel coupling [sJl- A de- 
tailed experimental study of these features can po- 
tentially provide a very interesting information on the 
details of strong interaction between charmed mesons, 
and might shed some new light on understanding 
'molecular'-type states. The phenomenon of reflec- 
tion of the threshold onset in one channel on another 
coupled channel is well known (see e.g. in the text- 
book |33|, Sect. 147). However the specifics of the 
pure P wave dynamics in the meson production am- 
plitudes in the e~^e~ annihilation and the apparent 
presence of resonance near the threshold bring in in- 
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teresting features, which justify a separate study of 
this phenomenon in the discussed context. 

In order to iUustrate the reflection of a resonance 
near the D*D* threshold we neglect here the isotopic 
mass splitting between the D* mesons and note that 
the resonance couples to a specific coherent mixture 
of the two possible production P-wave amplitudes for 
the D*D* channel: one with the total spin of the 
vector mesons S ~ Q and the other with S = 2. ^ 
Thus the final state of the D*D* pair can be consid- 
ered as one channel, referred here as H (heavy). A 
channel with a lighter meson pair, i.e. DD, D*D, or 
DgDs, is denoted here as L (light). The cm. mo- 
mentum in either of the L channels is already large 
and varies only slightly across the considered narrow 
energy range near the D*D* threshold, and this vari- 
ation is neglected here. 

Let W denote the total cm. energy E in the e+e^ 
annihilation relative to the threshold of production 
of the two vector mesons, each with the mass M: 
W — E — 2M, and let the 'nominal' position of the 
resonance be at the complex energy corresponding to 
W = Wq — ?ro/2. According to the general quantum- 
mechanical consideration, based on the unitarity and 
analyticity of the production and scattering ampli- 
tudes (see e.g. in the textbook .30|, Sects. 133 and 
145) at energy above the 'heavy' threshold, the reso- 
nant production amplitude Ah has the form 



W -Wo + ^iT^+glk^) ' 



(10) 



and the production amplitudes for each of the light 
channels are parametrized as 



Al = aL 



(11) 



where k = V MW is the cm. momentum of the D*D* 
pair. The coefficients aL,bL, and bH are complex pa- 
rameters, which are also taken to be constant in the 
present consideration. The magnitudes and the com- 
plex phases of these coefficients result from the rescat- 
tering among the 'light' channels and the correspond- 
ing absorptive parts associated with these channels. 
It is worth noting that there is no reason to expect 
the phase of the coefficient bn to vanish at A: = 0, due 
to the existence of 'light' inelastic channels strongly 
coupled to the 'heavy' meson pair. The normalization 
convention used here for the amplitudes corresponds 
to the cross section a{e^e~ H) — \Ah\'^ k for the 
'heavy' channel and cr(e^"e^ L) 



PL for each 



of the 'light' channels with being the cm. momen- 
tum in the corresponding 'light' channel at the 'heavy' 
threshold, which momentum is taken as a constant 
across the considered energy range. 

The expressions (I10|) and Hll() can be readily found 
from the standard summation of the 'blobs' in the 
graphs of the type shown in Figure 0] where the ul- 
traviolet part of the loop is absorbed into the defini- 
tion of the overall normalization and of the position 
Wq of the resonance. However the general treatment 
of the threshold behavior of the amplitudes guar- 
antees that these expressions are quite independent of 
assumptions about the form factors used in the graphs 
of Fig.H 



X 



X 



X 



L or H 



L or H 



L or H 



Figure 4: Tlie graplis for tlie self energy of the resonance 
X near the threshold for the 'heavy' channel (H). The 
absorptive and dispersive parts of the self energy arise 
from loops with the 'light' (L) and heavy meson pairs. 
The filled circle stands for the electromagnetic vertex. 

Furthermore, the 'nominal' width Fq is due to the 
coupling to the lighter channels, and the term k^ in 
the denominator in equations 1)10(1 and ((11(1 represents 
the additional contribution to the resonance width as- 
sociated with the decay into the 'heavy' meson pair, 
and gu is the coupling constant with the dimension 
of length. This extra contribution to the width can 
be written in terms of the excitation energy W as 
g'jjk^ = W yJW/e, where 



-3 



(12) 



is a parameter with dimension of energy. The critical 
point of the present consideration is that the ampli- 
tudes Al as functions of the energy W are analytical 
functions of W (with a cut at positive W) and the ex- 
pression ((11(1 can be analytically continued to negative 
values of W. Using the notation A = —W at < 0, 
the amplitudes Al below the 'heavy' threshold take 
the form 



A, 



aL 



(13) 



The term in the denominator in the latter formula, 
proportional to A'^/^, is clearly the analytical contin- 
uation of the g^ij k^ term in Eq. 1111 and the proper 
sign is as indicated in Eq. El'^- This term becomes 



^Another generally allowed amplitude, corresponding to S = '^This sign is due to the P wave motion, and is opposite 

2 and an F-wave, has to be small near the threshold and is to the one that would arise in a similar term for an S wave 
neglected here. (oc — A^/^). 
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more essential than the one hnear in A, starting from 
Awe, and gives rise to the discussed channel cou- 
pling effects. In particular, the expression ^Ay^— A 
cannot be arbitrarily negative and reaches its mini- 
mum, equal to -16e/27 at A = A™ = 16e/9. Thus 
in the case where the resonance is below the D*D* 
threshold, i.e. at negative Wo = — Aq, the real part 
of the resonance denominator has either two zeros if 
Aq < 16e/27, or no zero at all otherwise. In the lat- 
ter case the position of the minimum of the absolute 
value of the denominator, i.e. of the maximal effect 
on the cross section, is at A = A^, and does not de- 
pend on Aq at all, but rather is determined by the 
parameter e. The energy corresponding to A = A„i 
in fact also determines the position of the maximum of 
the resonance amplitude in the case where the 'nom- 
inal' position of the resonance is above the threshold, 
i.e. at positive Wq. This is due to the fact that at 
sufficiently small e the damping of the 'light' channel 
amplitude Al due to the absorptive term gjj in the 
denominator in Eq. ^2 suppresses the 'light' channel 
above the 'heavy' threshold. 

Clearly, the significance of the discussed threshold 
effect critically depends on the value of the parameter 
e for the near-threshold resonance. At least for two 
known resonances, V'(3770) and T(45'), at the thresh- 
olds of pairs of pseudoscalar heavy mesons the value 
of e can be found from the data: e « 60 MeV for 
?/'(3770), and e « 20 McV for T(45'). There is ev- 
ery reason to expect, from the observed onset of the 
D*D* production, that for the discussed possible res- 
onance the coupling to the vector meson pair should 
be stronger than for the case of pseudoscalar mesons, 
and the parameter e should be significantly smaller 
than for "0(3770). The parameter e is proportional to 
g]j'^ (cf. Eg. I12|l . Therefore, it is quite likely that for 
the D*D* threshold resonance this parameter can be 
smaller than that for the -0(3770) by a factor of 10 or 
more. It would come as no surprise if the value of e 
will be eventually found between 1 and 10 MeV. 

The possible types of behavior of the production 
cross section in a 'light' channel in the narrow energy 
range around the D*D* threshold is illustrated in Fig- 
ure|31for some representative values of the parameters, 
and where the isotopic splitting of the threshold is also 
taken into account. It is worth mentioning that the 
interference of the resonant amplitude with the non- 
resonant background is generally different for different 
channels, which may result in a significantly different 
behavior of the observed shape of the cross section 
for each channel. It is important to notice however 
that the features in the behavior of the cross section 
generally do not coincide with the position of the res- 
onance. The shape of the cross section of the e^e~ 
annihilation into the vector meson pairs D*D* above 
the threshold is only weakly sensitive to the param- 
eters of the resonance and is dominated by the very 



rapid rise from the threshold due to the P wave factor 



7. Angular Correlations in D*D* 
Production Near the Threshold 

As previously mentioned the threshold resonance 
should couple to a specific mixture of the states of the 
D*D* pair with the total spin S" = and S' = 2. The 
coefficients in this mixture and their relative phase 
can be studied by measuring angular correlations in 
the pro duction in e"*'e~ annihilation near the thresh- 
old [321 ■ It would also be of interest to study the 
behavior of these coefficients with increasing energy 
in order to possibly untangle the resonant and the 
non- resonant production amplitudes. 

Let ^0 s-nd A2 denote the production amplitudes 
for correspondingly the S = and S ^ 2 states, nor- 
malized in such a way that the total cross section is 
given by a{e+e- D* D*) = |^oP + |^2p. Then 
the distribution of the production rate in the angle 
between the direction of the momentum of either of 
the mesons and the e~'"e~ beams has the form 



da 
dcos( 



cx 1 



\A2? 



lOUol 



7|A2|2 + 10|Ao|2 



cos^ e . 



(14) 



Clearly, this angular distribution is not sensitive to the 
relative phase (p of the amplitudes Aq and A2. Such 
sensitivity however arises in the distribution over the 
angle between the momentum of the parent D* 
(D*) and the D (D) emerging from the decay D* 
Dtt {D* — > Dtt), or the similar correlation over the 
angle i)^ for the events with the decays D* 
[D* — > 1)7). The angular distributions are expressed 
through the two discussed amplitudes as 

(X |Aop -I- i- |A2p(13 + 21cos2 ^.) 
d cos Wtt 20 

4-^|^o||^2| cos(^(3 cos^t^^-l) , (15) 
v5 



and 



do- 



d cos 



cx |^oP + ^ 1^2^(13+ ;^sin2 ,9^) 



21 
Y 



-f— IA0IIA2I cos^(- sin^T^^-l) . (16) 

The latter two angular distributions however are not 
independent: one is found from the other by a sim- 
ple substitution cos^ ^ \ sin^i?^. The reason of 
course is that they both contain the same information 
about the vector meson polarization density matrix. 
Thus measuring one or the other, or both, depends on 
the specific experimental setup. 
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Figure 5: The types of behavior of the e'^e annihilation cross section (in arbitrary units) into one of the 'Ught' 
channels near the D* D* threshold. The horizontal axis in each plot spans the cm. energy range from 3980 MeV to 
4040 MoV. The D*"D*" and D*+D*' thresholds arc shown with shorter vertical tick marks. The plots are shown for 
three assumed values of the 'nominal' position Wo of the resonance X relative to the D'" D"'-'' threshold: -20, 3, and 20 
MeV, and the assumed value for this position is indicated in each plot by the longer vertical tick mark. In the plots in 
the upper row it is assumed that the coefficients Oi and hh are relatively real and have the same sign. In the lower row 
of plots the relative sign of these coefficients is assumed to be negative. The solid lines in the plots correspond to 
€ = 1 MeV, and the dashed one are for e = 10 MeV. The width parameter To is fixed at To = 50 MeV. 



8. Summary 

It is quite likely that in the resonance X(3872) we 
have come across a 'molecular' state made out of 
dominantly and D*^ . Although the existence of 
such states for "sufficiently heavy" heavy hadrons was 
not unexpected theoretically, it was not clear whether 
the charmed mesons were heavy enough. Given the 
present uncertainty in understanding the strong dy- 
namics of systems of heavy hadrons it is by far not 
clear yet whether the X(3872) is a result of coinci- 
dence of several factors facilitating the binding, and 
no other such states in the charmonium mass region 
exist, or there are other similar states with hidden 
charm. One such candidate for a molecular state is the 
apparent strong threshold singularity in e^e~ annihi- 
lation to D*D*. Further theoretical understanding 
of the strong interaction between heavy mesons can 
be (hopefully) facilitated by a more detailed study 
of the resonance X(3872) and by a detailed scan of 
the e+e" annihilation into exclusive channels with 
charmed mesons in the vicinity of the D* D* thresh- 
old. In this respect the interesting studies which 
look feasible include the measurement of the decays 
X(3872) -> D°D"tt, X(3872) D^D"j (the photon 
spectrum in this decay is directly related to the mo- 
tion of the mesons inside the X resonance), a search 
for the decay X(3872) D^D~^, which according 
to the presented discussion is sensitive to the tran- 
sition X(3872) ■!/'(3770)7. a search for the decays 
X(3872) XcjT^^ ■ For understanding the structure 
of the resonance very near the D*D* threshold in 
the e+e^ annihilation it is most interesting to have 
a detailed scan of the production cross section for the 



lighter pairs: DD, D*D and DgDs, since these chan- 
nels may exhibit a quite intricate behavior sensitive 
to detailed parameters of the resonance. Also a study 
of the angular correlations in e+e" D*D* near the 
threshold would provide an information on the 'hy- 
perfine' structure of the possible molecular state. 
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